Emittance fluctuation of mesoscopic conductors in the presence of disorders 
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We report the investigation of the dynamic conductance fluctuation of disordered mesoscopic 
conductors including ID, 2D and quantum dot systems. Our numerical results show that in the 
quasi-ballistic regime the average emittance is negative indicating the expected inductive-like be- 
havior. However, in the diffusive and localized regime, the average emittance is still negative. This 
disagrees qualitatively with the result obtained from the random matrix theory. Our analysis suggests 
that this counter-intuitive result is due to the appearance of non-diffusive elements in the system, 
the necklace states (or the precursor of necklace states in the diffusive regime) whose existence has 
been confirmed experimentally in an optical system. 
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The universal behavior of sample-to-sample fluctuation 
of mesoscopic conductors has attracted intensive theoret- 
ical and experimental studies in last two decades T]. For 
dc conductance it is well known that the quantum in- 
terference gives rise to a reproducible fluctuation in the 
coherent mesoscopic structures. When the sample size 
L is smaller than the coherence length but greater than 
the elastic mean free path, the system is in the diffu- 
sive regime and the conductance fluctuation exhibits a 
universal behavior with a universal conductance fluctua- 
tion (UCF) that is independent of Fermi energy, disorder 
strength, and system size. The numerical results of con- 
ductance fluctuation for various systems ^2] are in good 
agreement with the UCF values obtained from the dia- 
grammatic perturbation theory[l] and the random ma- 
trix theory 3]. For ac transport, due to the long range 
Coulomb interaction, the displacement current should be 
included in the ac current in a self-consistent manner i4!| . 
Hence it is important to examine the role played by the 
displacement current in the dynamic conductance fluctu- 
ation. For chaotic quantum dots, the fluctuating admit- 
tance and capacitance have been studied using random 
matrix theory [a, [Q]. Careful experiments have been car- 
ried out to measure this dynamic conductance in clean 
samples [7|, |8|. At low frequencies, the dynamic response 
of a system to the external bias is characterized by emit- 
tance which is the imaginary part of the low frequency 
admittance and is proportional to the partial density of 
state (DOS) of the systemjj]. Physically, for a conductor 
with large transmission coefhcient, the system responds 
like an inductor and the emittance is negative while for 
a conductor with low transmission coefficient the emit- 
tance is positive and the system behaves like a capaci- 
tor. For a ballistic conductor, the emittance E^ is found 
to be Ef, = -{\/4)e^dN/dE while for a metalhc diffu- 
sive wire the average emittance is Ef^ = {l/6)e'^dN/dE 
where dN/dE is the total DOS of the scattering region J9|. 
For a chaotic quantum dot, the emittance is also positive 
from the calculation of random matrix theory 5] showing 
capacitive-like response. For a disordered 2D waveguide. 



De Jesus et al have calculated the emittance distribu- 
tion using a continuous model [lOj. The emittance dis- 
tribution is Gaussian at weak disorder and crossover to 
a non-Gaussian distribution with a tail in the negative 
emittance region at strong disorder. 

For low dimensional systems, it has been predicted by 
Pendrvjll| that there exists quasi-extended states known 
as necklace states in strongly disordered regime. These 
necklace states are due to the multi-resonance and are 
very rare events. Recently, the existence of the analo- 
gous optical necklace states in disordered ID system has 
been confirmed experimentally in the Anderson localized 
regime|l2|. Since these necklace states come from multi- 
resonance, the electron dwells long time with large DOS 
in the scattering region giving rise to a large emittance. 
Thus we expect that the necklace states must have ma- 
jor influence on the average emittance for a disordered 
ID system. Motivated by this observation, in this paper 
we report an investigation of the disorder effects on the 
emittance and its distribution using tight-binding models 
for a carbon nanotube system (a quasi-one dimensional 
system), a two-dimensional waveguide, and a quantum 
dot. Our numerical results show that different from the 
previous physical expectation and theoretical results, the 
average emittance Ef^ is always negative even in the lo- 
calized regime. We attribute this counter-intuitive result 
as due to the non-diffusive elements in the system, the 
existence of necklace states (or the precursor of necklace 
states in the diffusive regime). 

Single wall carbon nantoube (SWNT)[13l| is a promis- 
ing candidate for the ac nanoelectronic transistor. In 
this work, we consider a disordered (6,6) SWNT of 480 
atoms connected to two semi- infinite perfect (6,6) SWNT 
leads. We employ the conventional nearest-neighbor 
tight-binding model whose Hamiltonian for the disor- 
dered SWNT reads H ^J2i e«4c^ + E,,jiU.j4cj + h.c), 
where c\ (ci) is the creation (annihilation) operator for an 
electron on the carbon site ?, e, = and i^j = SeT^ repre- 
sent the on-site energy and the nearest-neighbor hopping 
integral, respectively. Static Anderson-type disorder is 
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FIG. 1: The mean and rms values of emittance 
subband center E=3.2eV for different disorder strengths (one 
million configurations averaged). Inset shows the correspond- 
ing mean and rms values of conductance in unit 2e^ /h. 



added to e^ with a uniform distribution in the interval 
[—W/2, W/2] where W denotes the disorder strength. In 
this study, we have neglected the strong electron-electron 
interaction in SWNT. 

The conductance of the nanotube is calculated via 
the Landaucr-Buttikcr formula, G(0) = (2e^//i)T where 
T = Tr(riGT2C) is the transmission coefRcicent, G"" = 
[_E— i/— E^^ — Sj]"^ is the retarded Green's function of the 
central disordered SWNT and the relationship between 
self-energy and Hnewidth function is Fi^ = i[^i 2~^i 2]- 
With the help of density of states given by Green's func- 
tions, we can express emittance as [I,4j 



E,, = -Tr(Im[G"F2GTiG"], 
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where Da = CFaG", D^xx = [Da]xx, and [• • ■]^^ de- 
notes the diagonal element of the relevant square matrix. 
The quantities in the first and second terms in Eq.([T]) 
have clear physical meaning. The first term is the partial 
DOS dNi2/dE (up to a sign) describing the DOS of the 
electron coming from the second lead and exitingthe first 
lead while the second term consists of local D0S[9| and is 
due to the displacement current. The dynamic response 
given by E^ is either positive for small transmission co- 
efficient or negative for large transmission coefficient giv- 
ing rise to an inductive-like or capacitive-like behavior. 
For instance, for small transmission coefficent dNi2/dE 
is very small and the sign of the emittance is dominated 
by the second term which is positive. 

We first examine the emittance fluctuation defined 
by the root mean square (rms) values as rms(-E^) = 
li-^n) ~ i-^i^)^]^^^^ where (• ■ •) denotes averaging over an 
ensemble of samples, with different configurations of the 



same disorder strength. Fig.l shows the averaged emit- 
tance and its fluctuation vs disorder strength with the 
energy fixed at the subband center. We also plot the aver- 
age conductance (with unit 2e'^ /h) and its fluctuation for 
comparison. As the disorder increases the system crosses 
over from quasi-ballistic to diffusive regime. Eventually 
for strong disorders the nanotube becomes an insulator 
in the localized regime. The localization length £, defined 
as G = Gn exp(— 2L/^) is shown in Fig. 4a where Gn 
is the conductance in the absence of disorder and L is 
the system size. We see that the conductance fiuctuation 
in the diffusive regime W = [2, 4] (where ^ > L) ap- 
proaches to the one dimensional value of UCF 0.73e^/h 
(see the dotted line in the inset of Fig.1) obtained from 
the diagrammatic perturbation theory []| and the random 
matrix theory [3|. 

The main panel of Fig.l shows that the average emit- 
tance is always negative. This is a very counter-intuitive 
result since we expect that a conductor with small av- 
erage conductance gives a capacitive-like behavior corre- 
sponding to a positive emittance. It is true that for any 
specific configuration, the conductance and emittance are 
well correlated from our numerical result, i.e., large con- 
ductance corresponds to negative emittance and small 
conductance corresponds to positive emittance. The rea- 
son that the average emittance gives an inductive-like 
response is due to the so-called necklace states (or pre- 
cursor of necklace states in the diffusive regime) that are 
rare events in the disordered systems. Through these 
states, electrons can tunnel through the disordered sys- 
tem of size L via multiple scattering. The larger the 
system size, the smaller the number of necklace states 
is. The existence of these extended states in localized 
regime was predicted by Pendry|ll|. Different from the 
single resonant tunneling, the necklace state is a neck- 
lace of n quasi-extended states stretching from one site 
of the disordered sample to other. They occupy only a 
fraction of sites with fractal dimension 1/2 in ID rem- 
iniscent of percolation path[ll[. The analogous optical 
necklace states have been observed experimentally in dis- 
ordered multilayer svstems|12l| where the high transmis- 
sion peaks show up in the transmission spectrum while 
the average of logarithmic of transmission coefficient de- 
cays linearly with the thickness. Since the necklace states 
(multi-resonant states) are rare events, they have very 
long life time giving rise to a large negative emittance. 
To make this point clear, we plot in Fig. 2 the distribution 
functions P{Efj) versus emittance E^ at different disor- 
der strengths W . For this figure we have used data from 
200,000 different random configurations except in panel 
(h). We see that at low disorder strengths when W < 2 
the distribution is approximately Gaussian. In this dis- 
order range, as the disorder strength increases the values 
of emittance are all negative and the mean emittance 
gradually increases towards positive value. For larger 
disorder strength M^ = 4 in panel (d), the distribution 
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FIG. 2: The histograms of emittance at the subband center 
E=3.2 eV for different disorder strength W=0.5, 1.0, 2.0, 4.0, 
5.0, 6.0 and 10.0 eV, respectively. Panel (h) is the same as 
panel (g) but in a log scale of P(-E^) with one million config- 
urations. 



is still approximately Gaussian but the emittance can be 
positive while the mean emittance is negative. When 
W = 5,6 the distribution deviates gradually from the 
Gaussian distribution and the mean emittance is around 
zero but negative. At even stronger disorder W = 10, 
the distribution is sharply peaked at positive emittance 
with a long tail in the negative emittance (see Fig.2h). 
To be more precise, the distribution is asymmetric with 
large probability at small positive emittance while the 
probability of negative emittance is small but remains fi- 
nite even for large negative emittance resulting a negative 
averaged emittance. 

Why the effect of necklace states has not been seen in 
the conductance fluctuation but important for the dy- 
namic response? This is because at large disorders, there 
exists small number of necklace states with small proba- 
bility. For these necklace states the magnitude of emit- 
tance can be very large while the largest possible value 
of conductance is around one since it is unlikely for all 
multichannels to tunnel through. Therefore the contri- 
bution of necklace states to the average conductance is 
very small but is very large for the average emittance. 
For instance, in Fig.2h with W = 10, the average emit- 
tance is —0.11 and total conductance from eleven chan- 
nels is G = 0.24. Out of one million configurations, the 
total probability of these necklace states with emittance 
Ef^ < —100 is 0.3 percent. Their contribution to the aver- 
age emittance and average conductance are, respectively, 
E^ = -0.39 and G = 0.002. Glearly, these necklace 
states are important for average emittance and negligible 
for average conductance. Without these necklace states 
the average emittance would be positive. For exactly 



the same reason, the emittance fluctuation is greatly en- 
hanced due to the existence of these necklace states. We 
point out that these necklace states exist as long as the 
system size L is finite. In addition to the quasi-lD sys- 
tem, we have also studied the averaged emittance of a 
disordered ID tight-binding chain. Our result indicates 
that our conclusion remains, i.e., the average emittance 
is always negative for all disordered strengths. 

Note that the necklace states exist in the localized 
regime and our analysis above on these states is also done 
in the same regime. Now we wish to argue that the av- 
erage emittance E^ is negative in the diffusive regime 
as well. Since -E^ is a monotonic function of disorder 
strength, hence -E^ in the localized regime is larger than 
Efj_ in the diffusive regime. The fact that E"^ is nega- 
tive in the localized regime means that E^ is negative in 
all regimes including diffusive regime. Physically, in the 
diffusive regime, the system is roughly described by a dif- 
fusive conductor in parallel with multi-resonant channels 
connecting the leads (the precursor of necklace states). It 
is the appearance of these non-diffusive elements in the 
diffusive system, the precursor of necklace states that 
changes the sign of E^ in the diffusive regime from posi- 
tive to negative. 

Now we study the emittance using a tight-binding 
model for a quantum dot and a 2D system with dis- 
orders. The system size is E = 40a where a is the 
lattice spacing and two electrodes of width Ei are at- 
tached to the scattering region. For the quantum dot 
Li — E/4 while for the 2D system Li ~ L. The aver- 
age conductance and average emittance and their fluc- 
tuations for the quantum dot and 2D system are de- 
picted in Fig. 3 where we have collected one million con- 
figurations for each data point. From Fig. 3 we see that 
conductance fluctuation in the diffusive regime exhibits 
plateaus with universal values that are close to the the- 
oretical predictions UCF = 0.70e^/h for quantum dot'Sf 
and UCF — 0.86e^/h for 2D system^lj. For omittance, 
our results show that for both the quantum dot and the 
2D system, the emittances are negative independent of 
disorder strengths that disagrees with the theoretical pre- 
dictions. The emittance fluctuations for quantum dot 
(QD) and 2D system are consistent with ID case. 

Our results for the average emittance clearly disagree 
with the theoretical results obtained from the random 
matrix theory. This is because in order to use analytic 
approach, one has to make the following approximation 
on the second term in Eq.([T]) so that all the quantities 
involved are global. 



Trf 



DixxD2 



Dlxx + D2xx 



Tr(Di)Tr(D2) 
Tr(Di 



Da 



(2) 



It is easy to see that this approximation greatly over- 
estimates the second term in Eq.(IT]) at large disor- 
ders. For weak disorders with large transmission coef- 
ficient Dlxx and D2xx are in the same order of magni- 
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FIG. 3: The mean and rms values of emittance —E^ of the 2D 
system in panel (a) at _E = 0.62 and the QD in panel (c) at 
E = 3.0 for different W (one million configurations averaged). 
The corresponding mean and rms values of conductance in 
unit 2e /h are shown in panel (b) and (d), respectively, for 
the 2D system and the QD. 



tude hence Eq.([21) is a good approximation. However, 
at strong disorders there is a large mismatch between 
Dixx and D2xx for quite a number of sites x. Since 

DlxxD2xx / {Dlxx + D2xx) ~ D2xx whcn Dixx » D2xx 

the left hand side of Eq.([2]) is approximately given by 
Tr[min(Dixx, D2xx)] which is much smaller than the right 
hand side of Eq.Q at strong disorders. Using this ap- 
proximation, our numerical results are shown in Fig. 4. 
We see that, indeed, for both quasi-lD and 2D systems 
the average emittance with the approximation Eq.Q 
(solid square) overestimated the second term in Eq.(l) 
so that the results are very different from that obtained 
using Eq.(l). This is because the necklace states have 
been washed out in Eq.fl]). In addition. Fig. 4 shows that 
using approximation Eq.Q the emittance is positive in 
the diffusive regime in agreement with the theoretical 
predictions. 

In conclusion, we have studied the fluctuation of emit- 
tance and its distribution for disordered ID, 2D, and 
quantum dot systems. Our numerical results indicate 
that the average emittance is always negative indepen- 
dent of disorder strengths. Our findings disagree qual- 
itatively with the theoretical results obtained from the 
random matrix theory. The disagreement is due to the 
existence of non-diffusive elements, necklace states or the 
precursor of necklace states that are important for the dy- 
namic response and negligible for the conductance fluc- 
tuation. 

This work was financially supported by RGC grant 
(HKU 7044/ 04P) from the government SAR of Hong 
Kong and LuXin Energy Group. Endeavor Australia 



FIG. 4: The average emittance —E^j, and the localization 
length ^/L vs W (100 000 configurations averaged) for quasi- 
lD and 2D systems are shown, respectively, in panel (a) and 
(b). The other parameters are the same as in Fig.l and Fig. 3a. 
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